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Oi ■ A numerical and iterative approach for computing one-loop tensor integrals is 

presented. 

^H ■ 

> 

\£) ■ 1 Introduction 

o 

. _ . In Reference 1 a new approach has been introduced for computing, recursively 

f—) ' and numerically, one-loop tensor integrals. Here we describe a few modifica- 

tions of the original method that allow a more efficient numerical implemen- 
tation of the algorithm. We keep all of our notations as in Ref. 1 and, in par- 
ticular, put a bar over n-dimensional quantities and a tilde over e-dimensional 
O ,' objects (n = 4 + e). The formulae we want to modify are Eqs. (9), (35), (48), 

JL . (52) and (54) of Ref. 1 , that, all together, allow to reduce any (m + l)-point 

tensor integral with m > 2 to the standard set of scalar one-loop functions 2 . 
C^ , Such formulae are not symmetric when interchanging any pair of loop 

denominators D k , because are derived under the assumption that at least one 
of them (identified with Dq) carries a vanishing external momentum, namely 



H ■ D k = {q+p k ) 2 -ml, k = 0,---,m, p£ = 0. (1) 

Already after the first iteration, terms appear in which the denominator Dq 
is canceled by a Dq reconstructed in the numerator, so that the resulting 
integrals do not fulfill any longer the assumption of Eq. (1). A shift of the 
integration variable q is then needed to bring them back to a form suitable to 
apply the algorithm again. However, shifting q may generate a large amount 
of terms, especially when dealing with high rank tensors, so that deriving more 
symmetric formulae, in which p^ ^0, is clearly preferable. 
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A second useful modification is related to the problem outlined in Sec. 
6 of Ref. 1 , that occur when p\ = and p\ ^ (p\ ^ and p\ — 0) and 
(pi ■ V 2) ~ 0. For those kinematical configurations a new linear combination 
of the momenta p\ and p 2 is needed as a basis of the reduction to ensure the 
numerical stability of the algorithm. Once again, it is better to include such 
cases right from the beginning. 

2 The General Recursion Formula 

When^Q 7^ 0, the n-dimensional version of Eq. (9) of Ref. 1 should be modified 
as follows 

An) _ P_rp ( An)\a "I ,„. 
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where 

&% = (ftir 2 A + /i 2 ^)7^;.. T + (^+^)7^;p... T (l) (3) 

+ (r-i- + Cir 2 A ) 4"2i ;p ... T (2) - [r A (l + &) + r 2 A (l + &)] /<£* p ... r (0) . 

and where the extra integrals 7„^' p ... T are defined in Eq. (77) of Ref. 1 . 

In the previous Equations, fcj = j)j — po an d the massless 4-momenta ^2 
to be used as a basis of the reduction algorithm, as in Eq. (13) of Ref. l , are 
such that 

s 1 =£ 1 + ail 2 , s 2 = l-i + a 2 £i , (4) 

where si ]2 are suitable linear combinations of fci j2 

si = fci + £ifc 2 , s 2 = fc 2 + 6^1 • (5) 

By choosing 

6 = 2 si g n ( fc f) sign(fci • k 2 ) and 6 = ^ si S n ( fc i) si g n ( fc i ' fc a) , (6) 
the quantity 

7=- fifL_ s ( Sl . S2 )±VA, A = ( Sl - S2 ) 2 - S f S i, (7) 

(si -s 2 )t VA 



defined in Eq. (62) of Rcf. 1 only vanishes when k\ = k\ = {k\ ■ k 2 ) = 0, that 
always corresponds to collincar configurations cut away in physical observables, 
therefore solving the second problem outlined in the Introduction. The tensors 
T/xvXa , T^vx , 7j,„ and the 4- vectors r± 2 arc defined as in Ref . 1 , but in terms 
of ^i t 2 given in Eq. (4) and with the replacement p 3 — > k 3 . Finally 

h 2 = (m2-p£) + 6(m?-p?)-(l + $,)(rog-p§), 

^3 = (mj - p\) - {ml - pi) , 

- - ±-^- (8) 

The derivation of Eq. (2) closely follows the derivation of Eq. (9) of Ref. 1 . 
For example, choosing £i <2 as in Eq. (4), the quantity 

D li = ^[2(q-i 1 )i 2li + 2(q-i 2 )l llt ] (9) 

defined in Eq. (18) of Ref. 1 can be rewritten as 

D^ = [A + ^ 2 - (1 + ei)S + fti] r 2 ^ 

+ [5 2 + &£>i -(1 + &)A> + /» 2 ]ri M , (10) 

and generates the term j(") of Eq. (3). Analogously, choosing 6 = /c 3 in Eq. 
(22) of Ref. 1 , generates the first part of the last term of Eq. (2), because, 
when Pq ^ 

2(q-k 3 ) = D 3 -D Q + h 3 . (11) 

Finally, the equality 

q 2 = D Q + (ml-p 2 )-2(q-p )-q 2 , (12) 

is the origin of the second row of Eq. (2). 

3 Three-point Tensors 

When Pq =/= 0, rank 2 and rank 3 three-point tensor integrals need a separate 
treatment. The relevant formulae follow by adapting the theorems in Eqs. (37) 
and (40) of Ref. 1 to the case p% ^ 0: 
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The final results read as follows 
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4 Rank One Tensors 

In this section we adapt Eqs. (48), (52) and (54) of Ref. 1 to the case p^ ^ 0. 

4-1 The m — 2 case 

With t aa defined in Eq. (15) we get 
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^.,§ 27ie m — 3 case 

With T M „ A defined as in Eq. (23) of Ref. * we get 
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4-3 The m > 3 case 

The generalization of Eq. (54) of Rcf. l reads 
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hsW + 4^(3) - 4^(0) - ^ 3 a Ji n) A 



(18) 



where <5 = (4 • fc 4 )(4 ■ fc 3 ) - (£ 3 ■ k 3 )(£ 4 ■ k 4 ), and h 4 = (m| - p|) - (m§ - p 2 ,)- 
5 The Extra Integrals 

In most practical cases, the extra integrals, such as Im;p---r in Eq. (2), are 
either zero or scaleless, so that, even when p$ ^ 0, one can directly use the 
results given in Appendix B of Ref. 1 . In all other cases modifications are 
needed. For example, Eqs. (78) and (83) of Rcf. l must be replaced by 
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6 Conclusion 

We have derived a set of formulae to efficiently implement the n-dimensional 
reduction algorithm presented in Ref. . As for the three-point tensors, we 
limited our analysis to ranks < 3. For higher ranks, a general formula can be 
easily derived, with the help of Appendix C of Rcf. , using the theorems of 
Eq. (13). 
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